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FINAL  TECHNICAL  REPORT 

1.  Title:  Modeling  and  Parameterization  of  Gravity  Wave  Breaking,  Interactions, 

and  Instability  Processes  in  the  Lower  and  Middle  Atmosphere 

Agency:  Air  Force  Office  of  Scientific  Research 

Grant  Number:  F49620-92-J-0138  (also  AFOSR-91-0026) 

Principal  Investigator:  David  C.  Fritts 

Laboratory  for  Atmospheric  and  Space  Physics 
University  of  Colorado 
Boulder,  CO  80309-0392 

Program  Manager:  Lt.  Col.  James  G.  Stobie 

2.  Research  Objectives 

The  goals  of  this  research  were  to  address  two  sets  of  problems  relating  to  gravity 
w  ave  effects  in  the  lower  and  middle  atmosphere.  One  approach  w'as  intended  to  apply 
analytic  and  pseudo-spectral  collocation  techniques  for  studies  of  wave  excitation, 
instability,  dissipation,  and  interaction  processes  thought  to  play  important  roles 
in  the  circulation,  structure,  and  variability  throughout  the  atmosphere.  A  second 
effort  was  intended  to  apply  our  current  understanding  of  the  -spectfal  character  of 
the  atmospheric  motion  field  to  the  development  of  a  paraineterization  scheme  for  the 
most  important  wave  effects  in  large-scale  models  ofthe  atmosphere.  As  seen  below, 
these  goals  were  accomplished. 

3.  Research  Results 

The  contributions  made  in  each  of  the  major  areas  of  proposed  research  are  de¬ 
scribed  separately  below.  These  include  analytic  source  studies,  gravity  wave  insta¬ 
bility  modeling,  and  parameterization  development. 

a.  analytic  gravity  wave  source  studies 

Our  studies  of  gravity  wave  sources  used  Fourier  integral  and  Green’s  function 
techniques  to  assess  the  potential  for  gravity  w'ave  generation  by  sources  with  rel¬ 
atively  large  spatial  scales.  This  approach  was  taken  as  most  previous  efforts  have 
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focussed  on  excitation  of  gravity  -waves  at  relatively  smaller  scales,  despite  observa¬ 
tional  evidence  suggesting  a  clear  dominance  of  the  energy  within  the  gravity  wave 
spectrum  by  motions  with  large  scales  and  low  frequencies. 

Our  first  effort  addressed  excitation  by  geostrophic  adjustment  in  two  and  three 
dimensions.  In  each  case,  it  was  found  that  this  process  can  excite  waves  with  low 
frequencies  at  spatial  scales  in  reasonable  agreement  with  observations  near  the  s- 
tratospheric  jet,  where  much  of  this  activity  is  expected  to  occur.  These  results  were 
reported  by  Fritts  and  Luo  (1992)  and  Luo  and  Fritts  (1993). 

A  second  effort  addressed  a  source  that  has  been  the  subject  of  study  for  many 
years,  but  which  until  now  was  not  treated  in  a  quantitative  manner.  This  is  the 
possible  excitation  of  gravity  waves  by  a  solar  eclipse.  Our  study  defined  the  3-D 
wave  field  resulting  from  such  a  forcing  and  defined  the  altitudes  where  the  maximum 
responses  were  expected.  This  effort  was  described  by  Fritts  and  Luo  (1993). 

b.  modeling  achievements 

Under  AFOSR  support,  we  have  developed  2-D  and  3-D  nonlinear,  compressible 
pseudo-spectral  collocation  codes  in  association  with  Norwegian  colleagues  and  ap¬ 
plied  them  in  studies  of  gravity  wave  breaking,  instability,  and  turbulence  generation 
relevant  to  the  middle  atmosphere.  The  code  capabilities  and  the  results  obtained 
are  summarized  below. 

Our  initial  gravity  wave  breaking  simulations  focussed  on  a  wave  motion  with  a 
high  intrinsic  frequency  in  order  to  limit  the  computer  resources 'required  to  describe 
the  evolution.  We  also  performed  parallel  simulations  in  two  and  three  dimensions  in 
order  to  quantify  the  different  influences  on  the  incident  wave  amphtude  and  fluxes 
and  to  assess  the  validity  of  the  2-D  results.  Our  simulations  revealed  dramatic 
differences  in  the  two  evolutions,  suggesting  that  gravity  wave  instability  is  inherently 
3-D  and  that  2-D  models  are  unable  to  follow  the  physics  of  the  wave  breaking  process 
or  predict  the  influences  on  wave  amplitudes,  transports,  or  turbulent  mixing.  A 
comparison  of  the  2-D  and  3-D  results  and  an  assessment  of  the  instability  structure 
and  its  effects  were  provided  by  Andreassen  et  al.  (1993)  and  Fritts  et  al.  (1993). 
Vortex  breakdown  and  the  subsequent  evolution  toward  smaller,  more  isotropic  scales 
were  discussed  by  Isler  et  al.  (1993). 

Whereas  instability  in  the  2-D  simulation  was  confined  to  the  plane  of  wave  prop¬ 
agation  (artificially),  3-D  instability  w^as  found  to  comprise  counter-rotating  vortex 
structures  aligned  along  the  direction  of  wave  propagation  and  elongated  along  the 
imstable  phase  of  the  incident  wave  motion.  Relative  to  the  2-D  results,  the  3-D 
wave  structure  experiences  rapid  collapse  and  strong  amplitude  limits,  supporting 
the  convective  adjustment  hypothesis  and  suggesting  that  2-D  simulations  are  unable 
to  describe  the  physics  of  the  wave  breaking  process  or  the  imphcations  for  wave  and 
turbulent  transports.  Dramatic  differences  in  the  2-D  and  3-D  evolutions  are  also 
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apparent  m  the  transports  of  energy  and  momentum  and  in  the  evolution  of  eddy 
kmet.c  energy,  providing  further  evidence  that  2-0  simulations  are  not  suitable  for 
studies  of  gravity  wave  interactions,  instability,  or  induced  mixing  processes  in  the 
middle  atmosphere. 


c.  development  of  a  gravity  wave  parameterization  scheme 


An  important  component  of  the  research  effort  funded  by  AFOSR  was  the  devel¬ 
opment  of  a  spectral  parameterization  of  gravity  wave  energy  and  momentum  fluxes 
or  use  m  large-scale  models  of  the  lower  and  middle  atmosphere.  The  development 
o  his  scheme  and  its  applications  to  canonical  mean  wind  and  stability  profiles  and 

Lu  VanZandt  (1993),  Fritts  and 

Lu  (1993),  and  Lu  and  Fritts  (1993).  The  intent  in  this  development  effort  was  to 
create  a  parameterization  that  would  conform  as  closely  as  possible  to  observations 
of  the  atniospheric  gravity  wave  spectrum,  its  spectra]  shape,  and  its  variations  with 
dtitude.  This  we  beheved,  would  lead  to  more  realistic  fluxes  of  energy  and  momen¬ 
tum  with  altitude  than  predrcted  by  schemes  relying  on  discrete  and  noninteracting 

The  paraineterization  scheme  relies  on  a  spectral  description  of  the  atmospheric 
gravity  wave  field  that  is  largely  consistent  with  observations  throughout  the  atmo- 
sphere^  It  also  accounts  for  wave  filtering  and  anisotropy  in  a  manner  consistent 
with  theory  and  the  observations  that  are  currently  available.  As  a  result,  predicted 
energy  ^d  momentum  fluxes  are  c<^nstrained  realistically  arid  are 'in  broad  agree¬ 
ment  with  observed  mean  and  variable  values.  The  parameterization  was  designed 
e  implemented  efficiently  in  large-scale  models  of  the  middle  atmosphere  fit  is 

S  ACn  TIME-GCM  developed  by  Dr.  Ray  Roble.  in  the 

mod.^  P  f  atmosphere  chemistry  model,  and  in  several  middle  atmosphere 
_  els  in  Europe).  It  was  also  designed  to  incorporate  variable  wave  energies  and 
anisotropies  in  response  to  variable  gravity  wave  sources  and  source  strengths  in  the 
Wer  atmosphere,  but  this  variability  has  yet  to  be  determined  in  a  quantitative 

The  spectral  parameterizutiou  at  present  specifies  mean  gravity  wave  energies  at 
a  model  lower  boundary  and  computes  the  variations  in  wave  energy,  anisotropv, 
and  ener^  and  momentum  fluxes  with  altitude  as  functions  of  the  mean  wind  aL 
stabihty.  Tins  assumes  a  uniform  distribution  of  wave  sources  contributing  to  spectral 
energies  at  lower  levels.  Future  improvements  will  include  a  means  of  describing 

vSab  rr''"  77"^  specification  of  seasonal  wav: 

to  tidalfccbg  '“''S'  “'1  niomentum  fluxes  in  response 
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